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In this paper we study the influence of Dzyaloshinskii-Moriya (DM) interaction on quantum 
correlations in two qubit Werner states and maximally entangled mixed states (MEMS). We consider 
our system as a closed system of a qubits pair and one auxiliary qubit which interact with any 
one of the qubit of the pair through DM interaction. We show that DM interaction, taken along 
any direction (x or y or z), does not affect two qubit Werner states. On the other hand the 
MEMS are affected by x and z components of DM interaction and remain unaffected by the y 
component. Further, we find that the state (i.e probability amplitude) of auxiliary qubit do not 
affect the quantum correlations in both the states, only DM interaction strength influences the 
quantum correlations. So one can avoid the intention to prepare the specific state of auxiliary qubit 
to manipulate the quantum correlations in both the states. We mention here that avoiding the 
preparation of state can contribute to cost reduction in quantum information processing. We also 
observe the phenomenon of entanglement sudden death in the present study. 


I. INTRODUCTION 

Entanglement HEj and quantum discord [3] are the two aspects of quantum correlations and have the great potential 
in quantum information processing, which have many applications for future quantum technologies. The prominent 
role of entanglement have been investigated in teleportation, quantum secure cryptography, quantum games, quantum 
imaging, superdense coding and so on @H8j . The quantum correlations are very much fragile to the environmental 
interactions and they may loose their correlations. As a result they may lead to entanglement sudden death (ESD) 
for finite time EE HD]. The phenomenon of ESD in quantum information is a subject of investigation for practical 
quantum information processing. 

The two qubit system is the simplest bipartite system in quantum information that can be prepared either in pure 
or mixed bipartite quantum entangled states. Werner states are quantum states which are mixed two qubit entangled 
states m- Munro et al. have investigated another class of maximally entangled mixed states (MEMS) [121, which 
have more entanglement than Werner states and these are experimentally verified also m ■ Two qubit states have 
been realized in a variety of physical systems. So it is interesting to study the evolution of quantum correlations in 
physical systems under various interactions. In this paper we study the influence of Dzyaloshinshkii-Moriya (DM) 
interaction EHSl on the dynamics of entanglement and quantum discord HU in two qubit mixed bipartite states (i.e. 
Werner states and MEMS). We consider a pair of qubits and one auxiliary qubit, which interact with any one of the 
qubits of the pair through DM interaction. DM interaction plays a major role in quantum information processing and 
there is a vast literature in which the various quantum spin chains embedded in optical lattices governed by different 
models (Heisenberg, Ising, frustrated etc.) [IH, H] have been studied under DM interaction. DM interaction plays 
an important role to entangle, disentangle and amplifying the entanglement in quantum systems. But sometimes it 
may be the killer of entanglement in some quantum states. This has been shown by Zang et al. by considering two 
qubits A and B prepared in Bell states p20H!Z!2] and a third qubit C which interact via DM interaction with the qubit 
B. They have shown that DM interaction amplify and periodically kills the entanglement in two qubit Bell states. 
The entanglement manipulation can be done by controlling the state (i.e. Probability amplitude) of the third qubit 
(C) and DM interaction strength. Our observation tells that DM interaction is always not the killer of entanglement 
in quantum states as some of these remain unaffected by the same |23j . Recently, we also have shown that the state 
of auxiliary qubit (i.e. Probability amplitude) do not influence the entanglement in two parameter qubit-qutrit states 
|24j . In the present paper we show that the state of auxiliary qubit (i.e probability amplitude) does not come into 
picture to manipulate the quantum correlations in Werner states and MEMS; only DM interaction strength influences 
the quantum correlations. So in this case one can avoid to prepare the specific state of auxiliary qubit to manipulate 
the entanglement in both Werner states and MEMS. Avoiding the state preparation can help in cost reduction in 
quantum information processing. Recently the quantum networks have been established under DM interaction (25t - 
ehi, so the present study may be useful in manipulating and controlling the quantum correlations in quantum networks 
under DM interaction. 
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The plan of the paper is as follows. In Sect. 2 we present the Hamiltonian of the system. In Sect. 3 we discuss the 
two qubit Werner states and MEMS. In Sect. 4 we study the reduce density matrix (RDM) and eigenvalue spectrum. 
We study the influence of DM interaction on quantum correlations in Werner states and MEMS in Sect 5. Further, in 
Sect 6 . we present our conclusion. We present the quantum discord in Appendix and mathematica code to calculate 
the quantum discord for X-structured density matrices. 


II. HAMILTONIAN OF THE SYSTEM 

We consider a qubit (A)-qubit (B) pair and an auxiliary qubit (C) which interact with the qubit (B) of the pair 
through DM interaction. The Hamiltonian of the system can be framed as 

H = H ab +H™*, (1) 

where H AB is the Hamiltonian of qubit (A) and qubit (B) and H\ Jg, is the interaction Hamiltonian of qubit (B) and 
qubit (C). Here we consider uncoupled qubit (A) and qubit (B), so H AB is zero. Now the Hamiltonian becomes 

H = H' B q = D.{cfh x ac), (2) 

where D is DM interaction between qubit (B) and qubit (C), ah and d~c are the Pauli vectors associated with qubit 
(B) and qubit (C) whose components are Pauli matrices. We assume that DM interaction exist along the x, y and 
z-directions. In this case the Hamiltonian can be simplified as 


H x = D x .(a B (g> Gq - a B <g> cr£), (3) 

Hy = D y .(a% <g> a£ - a§ 0 a§), (4) 

H z =D z .(a B ®ac~a B ®ac), (5) 

where a B , a B and a B are Pauli matrices for qubit (B) and a@, Oq and a^, are Pauli matrices of qubit (C). The above 
Hamiltonians are the matrices having 4x4 dimension and easy to diagonalize by using the method of eigendecompo- 
sition. The unitary time evolution operator is easily commuted as 

U(t) = e~ iHt , ( 6 ) 

which is also a 4 x 4 matrix. This operator has been used to obtain the time evolution of the system. 


III. TWO QUBIT WERNER STATES AND MEMS 

In this section we discuss the two qubit bipartite Werner states (IT) and MEMS [12]. The Werner states in 2 x 2 
dimensions are widely used in quantum information and are invariant under unitary operators. The two qubit MEMS 
have been investigated by Munro et ah, which are more entangled than Werner states [12] . The entanglement in these 
states can be calculated by using the concurrence. The concurrence C in a density matrix p is given by 

C(p) = max{ 0, Ai — A 2 — A 3 — A 4 }, (7) 

where A*, i = 1,2, 3,4, are the square roots of the eigenvalues in decreasing order of pp?, pf is the spin flip density 
matrix given as 


pf = ( a v 0 a v )p*(a v 0 a v ). ( 8 ) 

Here p* is the complex conjugate of the density matrix p and a v is the Pauli Y matrix. The matrix ppf is used to 
calculate the concurrence. Two qubit Werner states can be written as 

Pw = + (1 ~7)j- (9) 

where is the singlet state and I is the 4x4 identity matrix. The state is entangled with 7 > 1/3 and carry the 
concurrence as C( 7 ) = (37 — l)/2. The MEMS are given below as 
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FIG. 1: Entanglement plot in Werner states and in MEMS 
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compare the concurrences in both the states as depicted in Fig. 1. In this figure Cw represent the concurrence in 
Werner states and Cmems represents the concurrence in MEMS. The concurrence achieves negative values for Werner 
states with 7 < 1/3, so it is disentangle within the range 0 < 7 < 1/3, while on the other the hand MEMS carry more 
concurrence than Werner states and hence more entangled than Werner states. 


IV. REDUCE DENSITY MATRIX (RDM) AND EIGENVALUE SPECTRUM 


To begin with we consider the auxiliary qubit in pure state as 

h P) = co|0) + c 0 | 1) 


which satisfy the normalization condition 


= 1 . 


( 11 ) 


( 12 ) 


In order to calculate the RDM first we calculate the time evolution density matrix of the system considered {i.e Qubit 
(A), Qubit (B) and auxiliary Qubit (C)}. It is given as 

P (t) = u(t) P m\t) 1 ( 13 ) 


where U(t) is the unitary time evolution operator obtained in Eq. (6) and p(0) is the initial composite density matrix 
of the open system (formed by qubit A, qubit B and qubit C). First we consider that the two qubits are prepared 
in Werner states, so the initial density matrix of qubits (A) and qubit (B) pair has to be considered as given in Eq. 
(9). We obtain the reduced time evolution density matrix by taking the partial transpose operation over the auxiliary 
qubit C. We observe that the factor (c§ + cf) has been involved with every element of the matrix, so by applying 
the normalization condition given in Eq. (12), the probability amplitude vanishes from the RDM and do not play 
any role in reduce dynamics. This result remains same when DM interaction is assumed in any direction (x or y or 
z). We obtain the same result when RDM is calculated for the MEMS. To obtain the RDM for Werner states with 
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FIG. 2: Entanglement plot vs time D x = 0.0,0.2, 0.4 and 0.6 


DM interaction in x, y and z directions we use the Hamiltonians H x , H y and H z given in Eqs. (3)-(5). Further we 
calculate the eigenvalue spectrum of RDM obtained for Werner states by considering all the Hamiltonians in order to 
calculate the concurrence by using Eq. (7). We obtain the same eigenvalue spectrum in x, y and z directions. The 
eigenvalue spectrum is given as 


{ j(-i + 7 ),+ 7), j(-l + 7 ),\{1 + 37)}- (14) 

Observing the eigenvalue spectrum we conclude that the components of DM interaction strength, D x , D y and D z in 
x, y and z directions respectively have not been involved. Hence the Werner states remain unaffected by the x, y and 
z components of DM interaction. Next we consider another initial state of two qubits (i.e. MEMS). First we consider 
the case with 2/3 < 7 < 1 and later with 0 < 7 < 2/3. 


A. Case 1: 2/3 < 7 < 1 

With the condition 2/3 < 7 < 1 for MEMS given in Eq. (10), we obtain the RDM by considering the DM interaction 
in x, y and z directions and calculate the eigenvalue spectrum Aj, i = 1,2, 3,4, for the matrix (fJfr ) to calculate the 
concurrence by using Eq. (7), we obtain the eigenvalue spectrum in x direction as 

{0,0, 7 , (-1 + 7 ) sin [2D x t] 2 }. 


(15) 
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FIG. 3: Variation of quantum correlations, QC (entanglement 


The parameter D x is involved in the spectrum, so the state is affected by the x component of DM interaction. In y 
direction the eigenvalue spectrum i = 1 , 2 ,3,4, of matrix (ppf) is given as 

{0,0,0, 7 }. (16) 

There is no parameter D y involved in the above spectrum, so the state is not being affected by y component of DM 
interaction. Finally in z direction we obtain the RDM and we get it as a X-structured matrix given below 
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The eigenvalue spectrum A i, i = 

1 , 2,3 

,4, of 

the matrix {ppf) is 

obtained as 



{ 0 , 0 , 7 , (-1+ 7 )sin[ 4 D-t] 2 }. (18) 

Observing the eigenvalue spectrum we conclude that the state is affected by z component of DM interaction. At this 
point we also mention that the quantum discord is easy to calculate in the matrix, as it is of X-structure type. 


B. Case 2: 0 < 7 < 2/3 

Under this case we obtain the RDM by considering the DM interaction in x, y and z directions for MEMS. First we 
obtain the eigenvalue spectrum A*, i = 1,2, 3,4, of the matrix (ppl) when DM interaction is along the x direction as 
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FIG. 4: Entanglement vs. time plot 


given below 


l ^/(3 — ^/5 — 4cos[4D a: I] — 2 cos [dlA^t]) 

{ 2 ’ °’ eTf ’ 

a/(3 + a/ 5 — 4cos[4D^tI — 2 cos[4D x tl) 

-- JVT -->• (19) 

We observer that parameter D x is involved in the eigenvalue spectrum and hence the state is affected by the D x . Next 
we calculate the same when DM interaction is assumed along the y direction. We obtain the eigenvalue spectrum A^, 
£ = 1,2, 3,4, of the matrix (ppf) as 


{^,°, 0 }. ( 20 ) 

Here we observe that there is no parameter D y involved, so the state is not affected by the y component of DM 
interaction. Next we obtain RDM and its eigenvalue spectrum, when DM interaction is assumed along the z direction, 
we obtain the reduce density matrix as given below 
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FIG. 5: Variation Quantum correlations, QC (entanglement and quantum discord) for D z = 0.0, 0.2, 0.4 and 0.6. 


_ 0 | cos[ 2 D z t] 2 

P ~ 0 -|sin[4 D z t\ 

.5 0 

The above matrix is X-structured, so quantum discord can be easily calculated. We obtain the eigenvalue spectrum 
Xi, i = 1, 2, 3,4, of the matrix (ppf) as follows 

{r,; g> °j 3 sin[AD z t]}. (22) 

Observing the eigenvalue spectrum we conclude that the state is affected by z component of DM interaction. 
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V. INFLUENCE OF DM INTERACTION ON QUANTUM CORRELATIONS 

In this section we study the influence of D x , D y and D z on MEMS for Case 1 and Case 2. We recall that Werner states 
are not affected by any component of DM interaction. So it is not possible to manipulate the quantum correlations in 


























two qubit Werner states by DM interaction strength and probability amplitude of auxiliary qubit. In this study we are 
able to calculate the quantum discord for MEMS only when DM interaction is along the z direction, because in this 
case we obtain that the density matrix is X-structured. Here we reiterate that MEMS are not affected by y-component 
of DM interaction ( D y ). Throughout this paper we plot the quantum discord by red color and entanglement by blue 
color. 

First we study the quantum correlation for case 1 (i.e. 2/3 < 7 < 1) in MEMS for different values of D x and D z . 
For D x we calculate the concurrence (Eq. 7) by using the eigenvalue spectrum given by Eq. (15). We obtain the 
entanglement in the absence of x-component of DM interaction ( D x = 0) and the result is shown in Fig. 2(a). We 
observe that as the value of 7 increases the entanglement increases linearly and achieve the maximum value 1. Now 
we increase the value of D x and plot the graphs for D x = 0.2, 0.4 and 0.6 in Figs. 2(b), (c) and (d) respectively. We 
observe that D x excite the entanglement and fluctuate it. Initially at t = 0, the entanglement achieves the value 0.67 as 
seen in figures, but as time advances it is excited and makes oscillations. As the strength of D x increases the frequency 
of oscillations of entanglement increases. Next we plot the effect of D z on entanglement and quantum discord for the 
same case. For this purpose we use the X-structured RDM given in Eq. (17). We plot both entanglement and quantum 
discord with different values of D z . First we calculate the quantum correlations in the absence of z-component of DM 
interaction ( D z = 0). We plot this result in Fig. 3(a). We conclude that the entanglement achieves the maximum 
value as 1.0, while the quantum discord achieves less value (0.95) than entanglement. Further we plot the quantum 
correlations for higher values of D z = 0.2,0.4 and 0.6 in Figs. 3(b), (c) and (d) respectively. Initially at t = 0, the 
entanglement achieves the value 0.67 but quantum discord achieves the value as 0.55. As the time advances both 
entanglement and quantum discord oscillates and quantum discord achieves less value than entanglement. When the 
strength of DM interaction increases the frequency of oscillations of quantum discord and entanglement increases. 
With the above results discussed we conclude that the component D x and D z do not produce the ESD in MEMS. 

Next we study the effect of D x and D z for the Case 2 (i.e. 0 < 7 < 2/3). First we calculate the concurrence 
(Eq. 7) by using the eigenvalue spectrum given by Eq. (19) for D x . We plot the effect of D x on entanglement for 
D x = 0.0, 0.2, 0.4 and 0.6 in Figs. 4(a), (b), (c) and (d) respectively. In Fig. 4(a) we plot the entanglement for 
D x = 0.0, in this case the entanglement achieves maximum value as 0.34. Observing the Fig. 4(b) corresponding to 
D x = 0.2 we get the result that at t = 0 the entanglement achieves the value 0.34, but as time advances the amplitude 
of the entanglement oscillates during the time interval 3.5 < t < 4.5, there is a gap in entanglement evolution (known as 
ESD). Further at D x = 0.4 we observe that the oscillations of entanglement increases and the entanglement smoothly 
goes to zero for a very small interval of time. But at D x — 0.6, the ESD zone almost vanishes and entanglement tends 
to zero and rise again. Next we calculate the quantum correlations for D z by using the X-structured matrix given in 
Eq. (21) and plot in Figs. 5 (a), 5(b), (c) and (d) for D z = 0.0, 0.2, 0.4 and 0.6 respectively. In Fig. 5(a) we plot the 
quantum correlation in the absence of DM interaction ( D z = 0) and conclude that the entanglement achieves higher 
value than quantum discord. Further we increase the strength of DM interaction ( D z = 0.2) and plot the quantum 
correlations in Fig. 5(b). Here we observe that quantum correlations fluctuate. At (t = 0) the value of quantum 
discord is 0.12, while the entanglement value is 0.35. Here we observe that quantum discord is less than entanglement. 
Again we find the phenomenon of ESD initially during the time interval of 1.8 < t < 2.1, but corresponding to this 
time interval the quantum discord is present. As the time advances the oscillations of both entanglement and quantum 
discord increases. In Fig. 5(c), with D z = 0.4 we observe that the entanglement smoothly goes to zero and rise again. 
But in Fig. 5(d) at D z = 0.6, we observe that the phenomenon of ESD vanishes and entanglement tends to zero and 
rise again. Here we mention that in the present case the states go under ESD by DM interaction during the short 
interval of time. It has also been noted that the frequency of ESD in MEMS for D x is half in comparison to the 
component D z . So the DM interaction in z direction is more intense than x direction. 


VI. CONCLUSION 

In this paper we have investigated the influence of x, y and z components of DM interaction on quantum correlations 
of two qubits. The qubits are initially prepared in Werner states and MEMS. We find that two qubit Werner states 
are not being affected by any component of DM interaction. While on the other hand the MEMS are being affected 
by the x and z components of DM interaction and remain unaffected by y component of DM interaction. Here we also 
observe that the state of auxiliary qubit (i.e Probability amplitude) do not play any role in manipulating the quantum 
correlation in both the Werner and two qubit maximally entanled mixed states. So one can avoid the intention to 
prepare the state of auxiliary qubit for the same purpose. The avoiding state preparation can help in cost reduction 
in quantum information processing. We also observe the phenomenon ESD in MEMS in the case 2 (0 < 7 < 2/3), 
while there is no ESD with the case 1 (2/3 < 7 < 1). When 0 < 7 < 2/3 the frequency of ESD is half induced by 
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x component of DM interaction in comparison to z component. So the DM interaction is more intense in z direction 
in comparison to x direction. Further we found that the quantum discord achieves less values than entanglement in 
MEMS with the case 0 < 7 < 2/3. 


Appendix A: Quantum discord 

To make this paper self-contained we briefly discuss quantum discord here m- It is known that a bipartite quantum 
state can carry classical and non-classical correlations. A theoretic measure of information is the mutual information in 
information theory. The mutual information stored in two classical distributions, equipped with two random variables 
X and Y, assuming valuesja;,:} and {y{\ with probabilities {pi} and {pj}, is defined as 

I{X : Y) = H(X) + H(Y) — H(X, Y ), (Al) 

or 


I(X : Y) = H(X) - H(X\Y ), (A2) 

where H{X) = —Y,iPilog 2 Pi and H(Y) = —Tiiqilog 2 <li are the Shannon entropies associated with random variables X 
and Y respectively, H ( X, Y) is the joint entropy of X and Y and H (A|T) is the conditional entropy of X and Y. From 
classical point of view both the above Eqs. (23) and (24) are equivalent but in quantum mechanical scenario both 
the expressions are different. Let us consider a bipartite quantum system equipped with parties A and B described 
by a composite density matrix p^ B with marginal density matrices p A and p B . The quantum analogues of mutual 
information of Eqs. (23) and (24) is obtained by replacing the Shannon entropies by von Neumann entropies as 

I(p AB ) = S(p A ) + S(p B ) - S(p AB ) (A3) 

and 


I(p AB ) = S(p A )-S(p A \p B ), (A4) 

where S(p A ) = — Tr(p A log 2 p A ) and S(p B ) = — Tr(p B log 2 p B ) are the von Neumann entropies of the subsystem A and 
B respectively and S(p AB ) = Tr(p AB )log 2 p AB is the joint von Neumann entropy of the composite quantum system 
AB. The quantity given by Eq.(26) is obscure, because the conditional entropy S(p A \p B ) depends upon measurement 
performed on either one of the party A or B. So both the Eqs. (25) and (26) are not equivalent from quantum 
mechanical point of view. The von Neumann projective measurement performed on system A projects the system into 
a statistical ensemble {p B ,p A }, such that 


where 


„ab 


Pk = 


(■ A k ®I B )p AB (A k ®I B ) 

Pk 


Pk = Tr[(A k 0 I B )p AB (A k 0 I B )\ 
A k = UM k U\ 
M k = \k){k\,k = 0,1 


(A5) 


(A 6 ) 

(A7) 


It is known that any unitary can be written up to a constant phase as 


U = tl + y.a, 

with t £ R, y = ( 2 / 1 , 2/2 j 2/3) G R 3 , cr is the Pauli vector of qubit and 


(AS) 


t 2 + y\ + y\ + 2/3 = 1 . 
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Now Eq. (26) becomes, with respect to to the projective measurement A k , as 


I(p AB \A k ) := S( P B ) - S(p AB \A k ), (A9) 

where S{p AB \A k ) := Y,p k S(p B ), is the conditional quantum entropy based on the measurement A k . Classical correla¬ 
tions is defined as the maximization of I(p Ab \A k ) over all possible measurements A k and is given by 

C(P AB ) ■= sup {Ak} I(p AB \A k ). (A10) 

On simplifying we get 

C{p AB ) := S{p B ) - min S{p B \A k ). (All) 

MO 

Quantum discord is given by the difference of mutual quantum information I(p AB ) given in Eq. (25) and classical 
correlations C(p AB ) given in Eq. (33). So we obtained the expression of quantum discord as 

Q(p AB ) = I(p AB ) - C{p AB ) = S(p A ) - S(p AB ) + min S(p s |A fc ). 

MO 

The calculation of quantum discord is based on complex maximization procedure and obtaining an analytic expression 
of C(pab) for general states is not an easy task. Here we write our x-structured density matrix as follows 


P 


AB 


Pl\ o 0 Pi4 
0 P 22 P23 0 

0 P32 P33 0 

. P41 0 0 P44 . 


(A12) 


Here we use the method to calculate the quantities Q(p AB ), CC(p AB ) given by C. Z. Wang et al.[25[|, So we obtain 
Q(p AB ), CC(p AB ) as follows 


(A13) 
(A14) 

where 

CCj = H(pn + P22) — Dj , 

QDj = H(pn + P33) + Y,'^_ l \ k Log2\ k + Dj, 

and 


CC{p AB ) = max(CCi,CC 2 ), 
QD(p AB ) = min(QD 1 ,QD 2 ), 


Di = H(t), D 2 = — ’Ej —1 Pjjl o<?2 Pjj - H(pn + p 33 ), 

H(pu + P 22 ) = ~{pn + p 22 )log 2 {pn + P 22 ) - [1 - (pn + P 22 )]log 2 [l - (pn + P 22 )], 

_ 1+'V / [ 1 — 2(p 33 +P 44)] 2 -t-4(| P 14 I + | P23 1 ) 2 
T — 2 > 

where X k are the eigenvalues of p AB and p 33 , P 44 , pn and P 23 are the elements of x-structured density matrix given 
in Eq. (35). 
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Appendix B: Mathematica code 


Under this section we provide the mathematica code to calculate the quantum discord for X-structured density 
matrices. The code is given below. 

m = {{pil, 0, 0, Pu}, {0, P22, P23, 0}, {0, P32, P33) 0}, {plA, 0, 0, Paa}}', 
e=Eigenvalues [nr]; 

Ai = e[[l]]; A 2 = e[[ 2 ]]; A 3 = e[[3]]; A 4 = e[[4]]; 

: 

Module [{xl = x}, —xl Log2[xl]-(l — xl) Log2[l — xl]] 

Log2[x_]=If[x = 0,0,Log2[x]]; 

p=FullSimplify[(ConrplexExpand [Abs[m[[l]] [[4]]]+ComplexExpand [Abs[m[[2]] [[3]]] ) 2 ]; 
r = (1/2) * (1 + y/{(l - 2(m[[3]][[3]] + m[[4]][[4 ]])) 2 +4(p))); 

DDl=N[H[r]]; 

a=-m[[l]][[l]]Log2 [m[[l]]m[[l]]]; 
b=-m[[2]][[2]] Log2[m[[2]][[2]]]; 
c=-m[[3]][[3]] Log2 [m[[3]][[3]]]; 
d=-m[[4]][[4]] Log2 [m[[4]][[4]]]; 
f=-H[m[[l]][[l]]]+m[[3]][[3]]; 

DD2=a+b+c+d+f; 

Qll == H[m[[l]][[l]]+m[[3]][[3]]]+Ai Log 2 [Ai]+A 2 Log 2 [A 2 ]+A 3 Log 2 [As]+A 4 Log 2 [A 4 ]+DD 1 ; 
Q22=H[m[[l]][[l]]+m[[3]][[3]]]+Ai Log 2 [Ai]+A 2 Log 2 [A 2 ]+A 3 Log 2 [A 3 ]+A 4 Log 2 [A 4 ]+DD 2 ; 

QD=N[Min[Qll,Q22]]; 

Quit]] 
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